Effects of spin and exchange interaction on the Coulomb-blockade peak statistics in 

quantum dots 
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We derive a closed expression for the linear conductance through a quantum dot in the Coulomb- 
blockade regime in the presence of a constant exchange interaction. With this expression we calculate 
the temperature dependence of the conductance peak-height and peak-spacing statistics. Using a 
realistic value of the exchange interaction, we find significantly better agreement with experimental 
data as compared with the statistics obtained in the absence of an exchange interaction. 
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The conductance through a quantum dot that is 
weakly coupled to leads displays sharp peaks as an ap- 
plied gate voltage is varied. Each conductance peak 
describes the addition of one more electron into the 
dot, and between peaks the conductance is "blocked" by 
the Coulomb interaction. The statistics of both peak 
heights and peak spacings in dots for which the single- 
electron dynamics is chaotic have been intensively stud- 
ied in recent years §. Some of the experimental obser- 
vations, e.g. the peak-height distributions at low temper- 
ature H, H, ^, could be explained at least qualitatively 
by the constant-interaction (CI) model, in which the in- 
teraction is represented in the simple form of an elec- 
trostatic charging energy. Other measured observables, 
such as the peak-spacing distribution , have indicated 
that spin and residual interactions beyond the charging 
energy should be taken into account. A consistent the- 
oretical approach that provides quantitative agreement 
with both the measured peak-height and peak-spacing 
statistics is still lacking. 
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Recently, a universal Hamiltonian was derived 
for a dot with a large Thouless conductance ^ 
(N is the number of electrons). An important contri- 
bution to the interaction part of this Hamiltonian is a 
constant exchange interaction in addition to the usual 
charging-energy term. The remaining interaction terms 
are suppressed at large gr- Here we study the effect of the 
exchange interaction on the finite-temperature statistics 
of both peak heights and spacings. To this end, we derive 
a closed expression for the conductance in the presence 
of a constant exchange interaction (in the sequential- 
tunneling limit). This formula expresses the conduc- 
tance in terms of quantities that characterize spinless 
non-interacting electrons. We then calculate the finite- 
temperature peak-height and peak-spacing statistics and 
find them both to be sensitive to the exchange interac- 
tion. Using an RPA estimate of the exchange interaction 
for the samples studied experimentally in Refs. 1^, ^, we 
obtain very good agreement with the observed temper- 
ature dependence of the standard deviation of the peak 
spacing. We also explain most of the known discrepan- 



cies between the experimental peak-height statistics g] 
and the predictions of the CI model for kT < 0.6 A (A 
is the mean spacing between spin-degenerate levels). 

The universal Hamiltonian of a quantum dot in the 
limit (7t oo is given by ||, 0| 
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where ex are spin-degenerate single- particle levels {a — 
±1 labels the spin). The second term in Eq. (|l|), where 
C is the dot's capacitance and fi is the total-particle- 
number operator, accounts for the electrostatic energy 
of the dot. The third term, in which S is the total- 
spin operator, describes a constant exchange interac- 
tion with strength Js- The occupation-number operator 
fix — nx+ + fix- of any single-particle orbital A com- 
mutes with the total spin, [nA,S] = 0, and the Hamil- 
tonian H is invariant under spin rotations. Thus the 
eigenstates of H are characterized by their particle num- 
ber N , the configuration of orbital occupation numbers 
n ~ {nx\ {nx = 0, 1 or 2), the total spin S, and the 
spin projection Sz = M. We label the eigenstates as 
\NTarySM) where the quantum number 7 distinguishes 
between states with the same total spin S and par- 
ticle configuration n. The eigenenergies are given by 

eis' = Ea CAriA + e'N^/2C - JsS{S + 1). 

In the limit of sequential tunneling (when a typical 
tunneling width is small compared with kT and A), the 
conductance can be calculated using a rate-equations ap- 
proach. In Ref. 1^, we developed such an approach in 
the presence of interactions and spin. In particular, an 
explicit solution exists when the orbital occupation num- 
bers nx are good quantum numbers. Expressing the con- 



ductance G in a rescaled form G 



^T/8hkT)g (where 



r is an average width of a level), we have, in the vicinity 
of the A^+l-st Coulomb-blockade peak 
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r^j^) are the single-particle 



Here gx = 2r 
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level conductances, where are the partial widths of 
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an electron in orbital A to decay to the left or right 
lead. The equilibrium probability of the dot to be in 
the state {Nn-ySM) is P^^^ = c-^^'^-s'-^fA')/^^ ^j^gj.g 
the partition function Z is a Boltzmann-weighted sum 
over all possible N- and {N + l)-body states (no other 
particle numbers contribute because of the charging en- 
ergy), and ep = + is an effective Fermi energy 
(ep is the Fermi energy in the leads, Vg is the gate 
voltage and ( = Cg/C with Cg the dot-gate capaci- 
tance). The Fermi-Dirac function f{x) = (1 + e^^)^^ 
is evaluated at an electron energy (relative to the Fermi 
energy) eg,g = e^,g, — e^^g — ep that conserves en- 
ergy at the transition between states |A^n7S'Af) and 
I (TV -I- l)n' j'S'M'). The corresponding reduced matrix 
element ([A^-|- l]n'7'S"||at^||Afn7S') enforces the selection 
rule S" =^ |S'± 1/2|. 

Eq. (^) can be rewritten in the form 

5 = E(^i°^+^i'^)5A. (3) 

A 

where the contributions with nx = and n\ = 1 are 
collected in w'f^ and respectively. For the cases 

with n\ — 0, the final (iV -f l)-particle state is given 
hy\{m)n'yS'M')^j:,,^{SM^a\S'M')aljNnjSM), 
where {SM^alS'M') is a Clebsch- Gordon coefficient. 
When n\ = 1 (and hence n\ — 2), the A^-particle state 
can be similarly related to the {N + l)-particle state by 
changing to a hole representation. This leads to the fol- 
lowing reduced matrix elements, 

(„Ys'i«ii„.5,.H--*{;^ ,4, 

Using the relation Piffi^'s) = Pi^s^'\^~I{4's)l we 
get 

=^T.^>^^^^sPn^s E (2^' + 1)/(4'5)> (5a) 

S S'=S±l/2 
I 



S' S=S'±l/2 

(5b) 

where the quantities b\,N,s = ^(("'A — l)(fT-A — 2))Ar^5 
and c\^N,s — ^{n\{n\ — l))iv,s ensure that the sum 
is only over contributions with = or 1, respec- 
tively. They are defined in terms of thermal expecta- 
tion values at constant particle number N and spin S, 

i.e. {X)n.,s = TrN,s[Xe-^"]/TrN,s[e'^"]- The quantity 
Pn.s is the probability to find the dot with N electrons 
and spin S, 

P^_s = e-''[^"'^+^"'^V^, (6) 

where Fn,s = — InTrAr^se^'^^^" '^"^■'"^'^ is the free 
energy of N non-interacting electrons with total spin S 
and Un,s = e^N^/2C - J,S{S + 1) - epN. 

The spin-projected trace of a scalar observable can be 
calculated from traces at fixed spin projection M using 
TrAr,5-'^ = ^rN^M=sX-TrN,M=s+iX. For spin-1/2 par- 
ticles, the projection on fixed particle number N and spin 
projection M is equivalent to projecting on a fixed num- 
ber of spin-up and spin-down particles n± — N/2 ± M. 
Therefore, 

Ttn.sX = tr«^_5 - tTN_^_g^^ N_^g^^^X , (7) 

where the traces "tr" on the r.h.s. are evaluated at fixed 
n+ and n_. Using X = e~^^ in Eq. (^ we 
find that the free energy in Eq. (^) is given by 

The free energy Fg in Eq. (||) is defined for q spinless 
particles e~^^'' — trqe^'^^^ ^ac^c^ where c\ and c\ create 
and annihilate spinless particles in non-degenerate levels 
with energies e\. The quantity c\^n,s from Eq. ( |5b| ) can 
now be expressed as 
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where fix is the particle-number operator of a non- 
degenerate orbital A. The function hx.N,s from Eq. (|5a| ) 
is expressed by replacing hx by [1 — fix) in Eq. (^. The 
complete expression for the conductance is then obtained 
from Eqs. (|), (|), (|), (|), (|) and the relation indi- 
cated in the previous sentence. Thus the dot's conduc- 
tance in model (|^) is determined in terms of the free 
energy Fq and single-particle occupation numbers {fix)q 



of q non-interacting spinless fermions. Both Fq and {fix)q 
are familiar from earlier works in the framework of the 
CI model, and can be expressed in closed form using 
particle-number projection [see Eqs. (140) in Ref. 

In chaotic dots, the single-particle Hamiltonian in (|l|) 
is described by random-matrix theory. We have stud- 
ied the statistics of peak heights and spacings for both 
the orthogonal and unitary symmetries. The dimension 
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FIG. 1: The width (t(A2) of the peak-spacing distribution for 
three different values of the exchange-interaction strength Js . 
The symbols are the experimental data of Ref. m\ . 



of the configuration sum in Eq. increases combina- 
torially with the number of single-particle orbitals and 
a direct use of (||) becomes impractical at higher tem- 
peratures. In contrast, the closed expression we derived 
greatly facilitates the calculation of the conductance for 
a rather large model system of 50 single-particle orbitals 
A. We checked that our results are not affected by the 
finite size of the system up to temperatures of kT '--^ 3 A. 

Theoretical calculations of the width a{A2) of the 
peak-spacing distribution, based on a spinless CI 
model describe qualitatively the observed decrease 
of this quantity with increasing temperature How- 
ever, a proper modeling of the peak-spacing distribution 
itself requires the inclusion of spin. When spin is included 
and in the absence of an exchange interaction, the calcu- 
lated values of cr(A2) (long-dashed line in Fig. |l|) show 
a large discrepancy with the experimental values (sym- 
bols). Fig. |l] also shows ct(A2) for non-zero values of 
Jg. For a gas constant of ~ 1.2 (that corresponds to 
the samples used in the experiments), the RPA estimate 
is Js ~ 0.3 A [|ll|, and we find for this value a very 
good agreement with the measurements. The results for 
Js = 0.5 A underestimate the experimental widths. We 
remark that at temperatures kT < 0.4 A, the model (|l|) 
does not describe well the shape of the peak-spacing dis- 
tribution, and it is necessary to include the fluctuating 
part of the universal Hamiltonian to explain the absence 
of bimodality At higher temperatures, the bi- 

modality is absent already in model (^ and the residual 
interaction has a negligible effect on the width. 

Another measured quantity is the ratio between the 
standard deviation cr((7niax) and the average ^max of the 
peak heights gmax [||- The experimental data for this 
ratio (symbols in Fig. ||) are seen to be suppressed in 
comparison with the results of model (|l|) without an ex- 
change term (long-dashed line in the left panel of Fig. |^). 
Spin-orbit interaction was proposed as a mechanism for 
this suppression at low temperatures p4[. It was neces- 
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FIG. 2: The ratio o"(gmax)/5max between the standard devi- 
ation and the average value of the peak height versus tem- 
perature kT. The left (right) panel shows data in the case 
of the elastic (rapid-thermalization) limit for three different 
strengths of the exchange interaction Js = (long-dashed), 
0.3 A (solid), and 0.5 A (short-dashed). The symbols are the 
experimental data of Ref. |8| . 



sary to assume a spin-orbit coupling that is sufficiently 
strong to completely decorrelate the spin-up and spin- 
down levels. However, spin-orbit effects are likely to be 
suppressed in the small dots used in the experiment. To 
determine whether an exchange interaction can explain 
the observed suppression of cr((7max)/3max, we calculated 
this ratio versus temperature kT for different strengths of 
the exchange interaction (see Fig. H). In the elastic limit 
(left panel), a realistic exchange interaction of Js = 0.3 A 
leads to closer agreement with the data. The remaining 
small discrepancy at temperatures kT < 0.6 A can prob- 
ably be accounted for by adding a realistic weak spin- 
orbit interaction. It still remains to explain the discrep- 
ancy at higher temperatures, where inelastic scattering 
may play a role. The calculation of Ref. |l^ showed that 
the suppression of o'(gmax)/ffinax due to inelastic scatter- 
ing is small for Js = 0. In the right panel of Fig. || we 
show results for the rapid-thermalization limit of strong 
inelastic scattering in the presence of an exchange inter- 
action. While the agreement (for Jg = 0.3 A) is now 
better at low temperatures, we do not expect inelastic 
scattering to be important at these temperatures. At 
higher temperatures, the rapid-thermalization limit does 
not describe the data, and it would be interesting to de- 
termine the effect of an additional weak spin-orbit term. 

For kT ^ A and Js — 0, the peak- height distribution 
-P(5max) can be calculated analytically ||] and is shown 
for the unitary symmetry as a solid line in the left panel 
of Fig. H (compared to the case of spinless electrons, the 
peak heights are rescaled [| by 8(a/2 - 1)^ - 1.37). Also 
shown (histogram) is the peak-height distribution calcu- 
lated at kT = 0.01 A and Js = 0.5 A. No significant 
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FIG. 3: Peak-height distributions. Left panel: The analyt- 
ically known distribution P((?max) at kT <C A and Js = 
(smooth curve) is compared with the corresponding distribu- 
tion at kT — 0.01 A and Js = 0.5 A (histogram). Right panel: 
Experimental data from Ref. [^] (gray-shaded histogram) at 
kT = 0.1 A are compared with the calculated distributions at 
Js ~ (dashed histogram) and Js = 0.3 A (solid histogram). 
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FIG. 4: The weak-localization parameter a versus tempera- 
ture kT in the elastic and rapid-thermalization limit for three 
different values of the exchange-interaction strength Ja. 



effect due to exchange is observed except for a small en- 
hancement of the probability at small peak heights. 

At finite temperature, the exchange interaction has a 
stronger effect on the peak-height distribution. The right 
panel of Fig. ^ compares the histogram (gray shaded) of 
the experimental data for P{g -max 1 9 max) at kT — 0.1 A 
with the calculated histograms for the cases of no ex- 
change (Js = 0) and Js = 0.3 A. This latter realistic 
value of the exchange interaction explains the observed 
suppression of the probability at small peak heights. 

The weak-localization effect in the average peak height 
attracted recent attention both in experiment and the- 
ory ||ia, 17, its suppression at higher tempera- 



tures was suggested as a signature of inelastic scattering 
in the dot. The effect is quantified by the parameter 
a = l — (g^^f/ftnajf)- III the rapid-thermalization hmit, 
a decreases rapidly with increasing temperature from its 
value of 0.25 at kT <C A In contrast, if inelas- 

tic scattering is negligible, a was expected to be tem- 
perature independent. However, calculations for Jg = 
showed a slight suppression of the elastic a around kT ~ 
0.25 A ||l|, |l|l. This was understood by the fact that 
close lying levels and hence higher conductances are more 
likely for the orthogonal symmetry. The effect of the ex- 
change interaction on a is shown in Fig. ^. We find that 
the dip in a around kT ~ 0.25 A is flattened out and in 
the small-temperature limit a becomes larger than 0.25. 
While a is seen to be sensitive to the exchange interac- 



agreement with experiment than in the absence of an 
exchange interaction. 
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tion at low temperatures, the experimental uncertainties 
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of Ref. |16| are too large to observe this effect. In the 
rapid-thermalization limit, a is insensitive to Js- 

In conclusion, we have derived a closed expression for 
the conductance in the presence of spin and exchange in- 
teraction. Using this formula we studied the dependence 
of the peak height and spacing statistics on the exchange 
interaction and found a significantly better quantitative 
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